Changing growth conditions during surface growth 
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Motivated by a series of experiments that revealed a temperature dependence of the dynamic 
scaling regime of growing surfaces, we investigate theoretically how a nonequilibrium growth process 
reacts to a sudden change of system parameters. We discuss quenches between correlated regimes 
through exact expressions derived from the stochastic Edwards- Wilkinson equation with a variable 
diffusion constant. Our study reveals that a sudden change of the diffusion constant leads to 
remarkable changes in the surface roughness. Different dynamic regimes, characterized by a power- 
law or by an exponential relaxation, are identified, and a dynamic phase diagram is constructed. We 
conclude that growth processes provide one of the rare instances where quenches between correlated 
regimes yield a power-law relaxation. 

PACS numbers: 81.15.Aa,68.35.Md,64.60.Ht,05.70.Np 
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I. INTRODUCTION 

Because of their ubiquity in nature, nonequilibrium 
growth processes have been the subject of numerous stud- 
ies during the last two decades P, 0, Q. Remarkably, 
these studies revealed very general properties of grow- 
ing interfaces that are encountered in a large variety of 
growth processes, ranging from crystal growth to tumor 
growth. In this context, due to its obvious technological 
relevance, thin film growth has been one of the major 
research thrusts. 

From the theoretical point of view, many important 
insights into the universal aspects of nonequilibrium 
growth processes have been obtained through the study 
of stochastic differential equations and of simple model 
systems [3,1- One of the simplest approach is due to 
Edwards and Wilkinson who described the surface 
growth due to particle sedimentation by dh{x,t)/dt — 
vV^hix, t) + u, where /i(x, t) is the surface height at time 
t at & site X of a d-dimensional substrate (of area .4) and 
u represents a constant flux of deposition. The physi- 
cal origin of the "diffusion constant" ly can be traced to 
the surface tension as well as T, the temperature of the 
substrate. When noise is added to this and the process 
is viewed from a co-moving frame of the steady state 
(i.e., h{x,t) ~ ut h{x,t)), we arrive at the stochastic 
Edwards- Wilkinson (EW) equation 



dt 



(1) 



where 7y(x, t) is a Gaussian white noise with zero mean 
and covariance (r?(x, t)r]{y, s)) = D6'^{x — y)S{t — s). A 
microscopic realization of the Edwards- Wilkinson equa- 
tion was soon proposed by Family d, 0| (see also Ref. 
Q for a recent more in-depth comparison of that model 
with the EW equation). In the random deposition with 
surface relaxation (RDSR) process particles drop from 
randomly chosen sites over the surface. Instead of stick- 
ing to the surface at the point of impact, the particles are 
allowed to explore the nearest vicinity of that point and 



are finally incorporated into the surface at the neighbor- 
ing site with lowest height. This diffusion step smoothes 
the surface and yields correlated growth. 

The roughness of a growing surface is characterized by 
the time dependent mean interface width 



W{t) = '^{{h-hf) 



(2) 



where h{t) = J h{x,t)d'^x/ A is the mean surface height 
at time t. In many instances the mean interface width 
(after an initial regime of random surface growth when 
starting from a flat initial condition) displays a power- 
law dependence on time, W{t) ~ , before saturating 
at a value W ^ where L is the size of the system. 
The growth exponent f3 and the roughness exponent a 
are universal quantities that characterize large classes of 
growth systems belonging to the same growth universal- 
ity class. Thus for the Edwards- Wilkinson universality 
class one finds for a one-dimensional substrate /3 = 1/4 
and a = 1/2. Other universality classes have also been 
found [1,07 some of which, as for example the Kardar- 
Parisi-Zhang (KPZ) and the conserved KPZ univer- 
sality classes [Ifl, , are of direct relevance for thin film 
growth. 

The morphology of growing structures can depend cru- 
cially on experimental parameters as for example the 
flux of deposited particles or the substrate tempera- 
ture p^ . Different experimental groups have reported 
a temperature dependence of the roughness of growing 
or sputtered surfaces, and temperature dependent values 
of the growth exponents have been found in some sys- 
tems [13l H m, m, mill m, m mmm- systems 

for which this has been observed include homoepitax- 
ial growth on Cu(OOl) [H, iH, Ag(lOO) 0, [1^, and 
Ag(lll) 16], amorphous thin- film growth '2^, growth of 
CdTe polycrystalline films [23], molecular-beam epitaxy 
growth of 81/81(11 11 jla] . as well as ion-sputtered 8i(lll) 
[i3, Ge(OOl) IWor Pt(lll) ^ surfaces. The de- 
pendence of the roughness on temperature can be rather 
involved, yielding different types of behavior for differ- 
ent systems, depending on how diffusion takes place and 
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whether additional smoothening mechanisms are present. 
In many instances one observes that growth becomes in- 
creasingly rougher with decreasing temperature, yielding 
an increase of the value of f3. On the other hand, some 
experiments [1^ [2^ revealed an increase of the global 
roughness with temperature. 

The observation of a transition between different dy- 
namic scaling regimes due to a change of experimental 
conditions is very intriguing and raises the question how 
a growing surface evolves from one regime to another af- 
ter a sudden change of, for example, the substrate tem- 
perature. We are not aware of any experimental study 
of growth processes where this kind of protocol has been 
implemented. However, we expect that the intriguing 
results reported here will motivate future experimental 
studies along similar lines. In this paper we discuss exact 
results derived from the EW equation with a variable dif- 
fusion constant, which allows us to investigate systemati- 
cally the changes in the surface roughness in case growth 
conditions are changed during the growth process. In- 
terestingly, different dynamic regimes are encountered, 
some of which are characterized by a power-law relax- 
ation. Here, we do not have a specific experimental sys- 
tem in mind. Instead, our interest is broader, namely, to 
understand universal responses of a growing surface to 
sudden changes of experimental parameters through the 
study of simple models. 

It is to be noted that our study is complementary to an 
earlier investigation due to Majaniemi, Ala-Nissila, and 
Krug '24]. Similarly to our work, these authors studied 
the impact of a change of growth conditions on processes 
described by linear growth equations. Whereas we dis- 
cuss in the following the effects of a variable diffusion 
constant, Majaniemi et al. analyzed how the growing 
surface reacts to a change of the noise in the system. 

There is also a second, theoretical motivation for our 
study. Sudden changes of external conditions, as for ex- 
ample due to a temperature quench, have been investi- 
gated in recent years in a large variety of systems, rang- 
ing from magnetic systems to glasses [1^, [2y, [23] . In the 
most common setting, a system initially prepared in some 
equilibrium state is suddenly brought out of equilibrium 
through a temperature quench. If the system is charac- 
terized by slow, non-exponential relaxation (as it is the 
case for a ferromagnet quenched to or below the critical 
point), interesting nonequilibrium processes and aging 
phenomena are observed |27[ . Some studies also focused 
on slow relaxation encountered in up-quenches in which a 
magnetic system initiall y in the ground state is quenched 
to the critical point [28ll29j. Interestingly, however, slow 
relaxation is usually not observed when both the initial 
and final temperatures are below the critical point. In- 
deed, for models with a discrete global symmetry such as 
Ising or Potts models a competing ordered state can not 
be reached if the starting point is close to one of the min- 
ima of the equilibrium free energy. Consequently, non- 
exponential relaxation is only encountered in this type 
of quench if the system has a continuous symmetry, as 



it is for example the case of the XY model [30l l3l| . As 
we will show in this paper, surface growth processes con- 
stitute an interesting class of systems where a change of 
parameters yields a transition from one correlated state 
to another characterized by a power-law relaxation. 

The paper is organized in the following way. In Sec- 
tion II we discuss the dependence of the surface width, 
derived from the EW equation, on the value of the diffu- 
sion constant. Section HI is devoted to the study of the 
time evolution of the surface roughness following a sud- 
den change of the diffusion constant. We thereby identify 
different dynamic regimes and present a dynamic phase 
diagram that summarizes the possible responses of the 
growing interface. Finally, we end with a summary and 
outlook in Section IV. 



II. INTERFACE WIDTH 

In order to study the effect of a change of external 
conditions on simple growth processes we consider the 
Edwards- Wilkinson equation ([1]) with a variable diffu- 
sion constant v (a microscopic realization of this process 
has recently been discussed in [13] )■ When describing a 
deposition or growth process by this equation, one implic- 
itly assumes that v and the noise amplitude D depend 
on the experimental parameters, as for example the tem- 
perature T. In the following we will not need to know 
the explicit dependence of v and D on these parameters 
(which would be system dependent), and study how the 
interface width changes when changing the value of v. 
The reaction of the growing surface to a change of the 
noise has been studied in [24| . 

The stochastic EW equation ([T]), can be solved exactly 
to give us for a fixed value of v the width (squared) 



1 - e 



2vL ^ 



(3) 



where fc„ = l-KnjL and the sum is over ]— L/2,L/2] but 
excluding the zero mode: n — Q. In Figure [T^,c we show 
the time dependence of the surface width for, respec- 
tively, the case with fixed L = 1000 at different i^'s and 
the case with a fixed v = 0.1768 and various L's. As for 
the RDSR process one distinguishes three regimes sepa- 
rated by two crossover points: a random deposition (RD) 
regime, followed by a EW regime, with a final crossover 
to the saturation regime. In contrast to Family's origi- 
nal model, the initial RD process is not confined to very 
early times t < \ but might extend to larger times. In 
fact, the crossover time t\ between the RD and the EW 
regimes is shifted to higher values for decreasing diffu- 
sion constants and diverges in the limit of vanishing v. 
As the crossover is smeared out, we identify the crossover 
point with the intersection point of the straight lines fit- 
ted to the two linear regimes in the log-log plots. We 
have the identity W"^ = t in the RD regime, yielding 
the width W\ = \/t\ at the crossover point. In the EW 



3 



regime the relation between width and deposition time 
changes toW cc t^/^. This regime extends up to a second 
crossover point {t2^W2), whose precise location depends 
on the values of v and D and beyond which the final 
saturation regime prevails. The crossover between the 
different regimes is further illustrated in Fig. WpA where 
we show the time evolution of the effective exponent 
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dlogW _ 

dlogt ~ ^ [1 _ /fc2 



(4) 
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FIG. 1: (a) Log-log plot of the surface width vs time for a sys- 
tem of size L = 1000 and different diffusion constants. The 
dotted lines have the slopes 1/2 and 1/4 expected in the ran- 
dom deposition and EW regimes, respectively. The locations 
of both crossover points depend on the diffusion constant. 
The data are obtained from the exact solution of the EW 
stochastic equation, (b) Time evolution of the effective ex- 
ponent ((U for the data shown in (a), (c) Log-log plot of the 
surface width vs time for systems of different sizes evolving at 
the same diffusion constant v = 0.1768. (d) Time evolution 
of the effective exponent Q for the data shown in (c). The 
data in this and the following figures have been obtained for 
D = l. 

We end this Section with a few remarks: 

• The ("first") crossover from the RD to the EW 
regime, denoted by {ti,Wi), depends only on i/. 
For the range of L's we explored, we find Wi = 
ti = r/h', with a constant r = 0.148. In the i — > oo 
limit, T ^ l/27r (as shown in the Appendix, see also 

(H). 

• The ( "second" ) crossover from the EW to the sat- 
uration regimes, denoted by (^2,^^2)1 depends on 
both ly and L. W2 may be identified with the sat- 
uration width, i.e., (^DL/8n^i')J2n^~^ t23i]. Com- 
bining this result with the line drawn through the 



EW regime, we arrive at ^2 
pendix) . 



(L/24r) ti (see Ap- 



III. SUDDEN CHANGE OF GROWTH 
CONDITIONS 



With a clear picture of the properties of a surface grow- 
ing under a constant diffusion constant, we proceed to 
study the time evolution of the roughness when v is sud- 
denly changed. To investigate the response of the growth 
process to this change, we use the following protocol: We 
start at i = with a flat surface and let the surface grow 
at Vi until time t = s, at which point we change the 
diffusion constant to the final value Vf. The change of 
roughness is then monitored through the time evolution 
of W^. 

As our system displays three different roughness 
regimes (RD, EW, and saturation), there can be in prin- 
ciple nine scenarios for the change to be arranged. They 
can be distinguished conveniently by (1) i>i, the diffusion 
constant of the initial growth, (2) Vf, the final value of 
the diffusion constant, and (3) s, the time at which 1/ 
is suddenly changed. By choosing these three controls 
judiciously, we can access all the scenarios. However, 
covering all cases in detail is not the aim of this paper. 
Instead, we are interested in the new phenomena asso- 
ciated with the ^-dependence of the width, W'^{t,s), af- 
ter (up- or down-) quenches. Obviously, for t ^ s, the 
surface roughness will settle into the value in an "unper- 
turbed" system, grown at from the start. Denoted by 
W^{t), it is just expression ([1]) with 1/ = i/f. We also 
refer to this as the "reference system." To highlight the 
changes, we also study the difference (with t > s) 



AW^it,s) = \W^it,s)-Wlit)\ 



(5) 



between quenched system and the reference. Clearly, this 
quantity reveals how the roughness of the growing surface 
adapts itself to the new "experimental" condition, and 
behaves very differently for the various cases. Our goal is 
to map out the regions in the Ui-Vf-s space corresponding 
to the novel behavior following a quench. 

To be precise, we will evolve the height h{x, t) starting 
from h{x, 0) = with Vi to time s and then continue 
with Vf until time t. At that point, we compute the 
width squared and denote it by W'^{t, Vf] s, Vi). 

Our starting point is the exact solution of (P) 



h{K,t) = e~'''''r.(t-s) j dt'e---'=-.(--*')^(fc„,i') 


t 

+ /dt"e-'^/'="(*-*")77(fc„,t") (6) 



written here in terms of the Fourier amplitudes for h{x, t) 
and ri{x,t): h{kn,t) = J dx e'''^"^ h{x , t) , etc. Since 
the noise is delta-correlated, simplifies so that the 
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width square is 

W^{t,lyf■,S,I^^) ■ 



D 

T 



E 







(7) 



Since the width square W^{t, vf) of the unperturbed sys- 
tem is given by 



D 

T 



z 



-2!/j:feJ(t-t') 



(8) 



we arrive at the exact resuh for the difference AW'^ 
\W^it,i.f;s,,.,)-Wi{t,,.j)\ : 



D 



-2i^/fc„(t-s) 



dt' 



D 



E' 



'2vfk1^(t-s) 



-2uikts 



(9) 



For later convenience, we define = VfAW'^{t, Vf \ s, Vi) 
and note that it depends only on three scaling variables: 



tJL=v.Jvf, a = VfS, p = tls . 
Explicitly, we have 



E' 



-2fc^cr(p-l) 



-2k^fj,iT 1 — e~ 



2fe., 



2kl 



(10) 



(11) 



To set the stage for discussions, we begin with the data 
for some typical cases, all with s = 10^, shown in Fig. 
[5J In order to be able to discuss the different cases for 
a fixed s, we must work with a relatively small system: 
L = 400. The dashed lines in (a-c) represent W^{t) in 
an unperturbed system. With i^f = 4.7 • 10~^°, 0.00047, 
and 0.18, the surface is, at the time of the quench, in the 
(a) RD, (b) EW, and (c) saturation regimes, respectively. 
The corresponding differences, AW'^{t, s), are shown in 
Fig. Hi-f. 

The effects of two up-quenches into the RD regime, 
from the EW (m, = 0.0046) and the saturation {vi = 0.23) 



regimes, are displayed in Fig. [2^. As Fig. [2ji shows, for 
quenches into RD, the width W^{t, s) cannot reach that 
of the reference system, Wi{t). The "best" AW^{t,s) 
can achieve is a constant. The physical origin of this be- 
havior lies in the linear growth of in the RD regime. 
Thus, for two unperturbed systems started at different 
times (say, t = to and ti), the difference Wq — is just 
a constant: ti — to. In our case, the correlated growth up 
to time s endowed our surface with a smaller W"^ (s) than 
the reference W^{s). Immediately after the quench, cor- 
related growth is simply replaced by independent growth 
of different columns and the width (s) is "frozen" in 
as a kind of "initial condition" (at t — s). As a result, 
the difference AW^ = W^{t) — W^{t,s) remains at the 
value W^{s) — (s). Of course, if we follow these two 
systems further in time, AW'^ will eventually vanish. 

Turning next to quenches to the EW regime, we found 
the most interesting behavior (Fig. [2l3,e). The difference 
AW^ initially decreases rapidly, before crossing over to 
a slower, power-law decay at larger times (t 3> s): 



AW^ 



(12) 



For example, for the cases shown in Fig. [2^, we mea- 
sure close to the end of the time interval the exponents 
7 = 0.72 for ly, = 0.23, 7 = 1.24 for = 0.0046, and 
7 = 1.65 for — 4.7 • 10^^°. As we argue below, these 
are effective values, the asymptotic values of 7 being 1/2 
or 3/2. Below we will also discuss in more detail the 
conditions under which these values can be expected. Fi- 
nally, for a quench to the saturation regime, AW"^ decays 
exponentially: 



AW^ - exp(-Kt) 



(13) 



with a decay constant k that depends both on the value 
of the final diffusion constante Vf and on the system size 
L. 

Up to now, we have shown only the simplest situation 
where the system is well within a given initial regime at 
the moment of the quench and, in addition, that it has 
time to relax into a well defined regime of the reference 
system. Clearly, as we let the final system evolve further, 
it may crossover to a different regime (e.g., in case of Fig. 
[2^,d a crossing over to the EW and saturation regimes 
will take place for larger t). Therefore, we should expect 
the general relaxation process to be quite complex. 

The closed forms (|9llip are not particularly transpar- 
ent, as they involve all possible crossover behaviors. To 
shed some light on the various scenarios, we consider 
some limiting cases where simple properties (exponen- 
tials and powers) can be extracted. A straightforward 
case is i'f{t — s) ^ L^, so that i^/fc^(t — s) is always 
large. Then, the leading decay is exponential, namely 
g-STT v]{t-s)/L ^ sij-^ce the other terms in the sum will be 
much smaller: 



ui{t~s)/L' 



(14) 
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FIG. 3: Contour plot of 7 as a function of UfS and i^i/i^f for 
t/s — p = 64. Four different regimes, separated by crossover 
regions, are identified. The two dashed lines separate the 
three qualitatively different types of behavior encountered 
when plotting the effective exponent as a function of t, see 
Fig. [S] 



FIG. 2; (Color online) (a-c) Time evolution of the width 
square in case the diffusion constant is changed after 10^ time 
steps. The dashed lines show W^it) for an unperturbed sur- 
face growing at constant i/f. In (a) the quench is to the RD 
regime, whereas in (b) and (c) the quenches are to the EW 
and saturation regimes, respectively, (d-f) The same cases as 
shown in (a-c), but now the difference AVK^, see Eq. ([5]), is 
plotted. Qualitative different behavior is observed, depending 
on the regime that the unperturbed system has at the quench 
time. 



In the opposite limit, where i^fS, ViS <^ 1, we have 
"1 - e^^"'''"* l-e'"^""!^" 



2vfkl 



^{yf-v,)s'ki (15) 



to leading order. The summation over n then yields, for 
t — 1, a, power-law decay: i.e., 7 = 3/2. 

In order to explore the three-dimensional parameter 
space (/i, (T, p) in a comprehensive way, we evaluate nu- 
merically the closed form (jlip . The exponent 7 can be 
defined effectively as —d\og{ri) / d\og{t). In Fig. [5] we 
show the contour plot of 7 as a function of /i and a for 
p — 64. This plot reveals four different regimes: the 
regime where Q, or is constant (labeled by 7 = 0), 

two power-law regimes with values 7 = 1/2 and 7 = 3/2, 
and finally a regime of exponential decay for large cr. 
The different regimes are separated by crossover regions 
where the effective exponent does not lock-in into one of 



the values 0, 1/2, or 3/2. 

Fig. 11^ shows how the extensions of the four regimes 
depend on the value of p in cases where p ^ 1. In- 
terestingly, an increase of the value of p mainly shifts 
the contours in the \og{(j) vs. \og{p) plot along the 
(—1, 1) direction. This is shown in Fig. where we plot 
\og{v ft) = log(crp) vs. \Qg(v^s/vft) = \og{p/p). This 
way of plotting indeed leads to an approximate data col- 
lapse, which gets better for larger values of p. That this 
data collapse is only approximate also follows from in- 
spection of the exact solution (fTTjl . Still, Fig. |4|3 nicely 
allows us to visualize the extent of the different dynamic 
regimes for large ratios p. 

Finally, in Fig. [5l we discuss the change of the effective 
exponent 7 as a function of t for various values of p, (and 
s = 10^). Note that an increasing time t corresponds in 
Fig. [3] approximately to a cut along the (—1, 1) direction, 
so that we can distinguish three typical scenarios, sepa- 
rated by the dashed lines there. Along the upper dashed 
line. Fig. [5^ shows the effective exponent rising to the 
7 = 1/2 plateau where it remains for a long time before 
crossing over to the regime where the difference AW^ 
vanishes exponentially fast ("7 = 00"). For the region 
above this line in Fig. [31 we can expect similar results. 
Along the lower dashed line, the same behavior is seen, 
except that the plateau value is now 7 = 3/2 (Fig. [St). 
This can also be expected for the region below this lower 
dashed line in Fig. [3) Between these two protocols, a 
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FIG. 4: (Color online) (a) Contour plots of 7 as a function 
of i/fS = a and Vi/i'f for t/s = p = 64 (full lines) and 256 
(dashed lines). Only contours bounding the 7 = 0, 1/2 and 
3/2 regimes are shown, (b) The same contour plots as shown 
in (a) but as a function of Vft and Visjvjt. An approximate 
collapse of the contours is observed. 

more complex behavior is encountered, as the effective 
exponent shows some tendency to lock-in at both val- 
ues 1/2 and 3/2 (see Fig. [5)d). We should remind the 
reader that the 7 (t) curves shown here are applicable for 
all varieties of quenches (different quench times, as well 
as different initial and final values of v) as long as the 
rescaled time i^/s is fixed. 

It is interesting to note that a change of the noise dur- 
ing the growth process also yields a power-law relaxation, 
as discussed in [23|. Thus for the one-dimensional EW 
equation it was observed that AVF^ ^ t~^/^ if the system 
is in the EW regime before and after changing the noise 
of the system. The situation studied in 241 is therefore 
comparable to our case, even though we do observe a 
richer behavior when changing the diffusion constant. 



IV. DISCUSSION AND OUTLOOK 

The morphology and roughness of growing surfaces de- 
pend in a crucial way on experimental conditions. Moti- 
vated by various observations of a change of the rough- 
ness universality class when changing experimental con- 
ditions, we have studied in this work how a growth pro- 
cess reacts to a sudden change of the diffusion constant. 
Exploiting the fact that the stochastic EW equation can 
be solved exactly, we carried out a comprehensive study 
of the response of the growing interface. Four main re- 
laxation regimes, separated by crossovers, were identified. 
For extremely long times any finite system will eventu- 
ally relax in an exponential way, due to the presence of 
the saturation regime, but this only takes place after an 
earlier power-law relaxation. For finite times, a growing 
surface rapidly reacts to a change of v such that its mor- 
phology (roughness) approaches that of a reference sys- 
tem that was allowed to grow under stable external con- 
ditions. If the change of growth conditions takes place 
at a time where both the quenched and the reference 
system are in a correlated growth regime, than the relax- 
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FIG. 5: The three different quench types illustrated by the 
time dependence of the effective exponent 7. In all cases the 
quench takes place at s = 10^. Type (a) includes quenches 
with initial conditions VjS and Vi/vf located above the upper 
dashed line in Fig|4l For type (b) the starting point is located 
between the two dashed lines in Fig 3] wheres for type (c) 
the initial conditions put the starting point below the lower 
dashed line. 



ation process is governed by power-laws over many time 
decades. This is our main result (see also [24j]), and we 
expect this to hold true for other systems. 

Of course, we do not expect that our results for a 
one-dimensional system can be applied quantitatively to 
any of the experimental systems in which a tempera- 
ture dependence of the roughness exponent has been ob- 
served. Nevertheless, we expect that these results should 
be generic for growth processes with a sudden change of 
external conditions, and that the intriguing signatures 
revealed in our study can be generalized to more realistic 
models, so that they can be observed in experiments on 
physical systems. 

From a more theoretical point of view, our study 
reveals that growth processes provide one of the rare 
cases where a power-law relaxation is generically ob- 
served when quenching from one correlated state to an- 
other. More specifically, the exactly solvable Edwards- 
Wilkinson equation allows us to derive analytical expres- 
sions, thus permitting a complete investigation of the 
various possible scenarios. It is to be expected that the 
power-law relaxation encountered at a quench also en- 
tails interesting aging processes not studied in the past. 
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Indeed, in the few published studies of aging in growth 
processes d, HI, constant model parameters were al- 
ways assumed. 

Our study can be extended in various directions. On 
the one hand, we can study quenches in systems where 
the interface is stabilized not by surface tension, but by a 
curvature Hamiltonian. The simplest case is given by the 
noisy MuUins-Herring equation [10, 34,]. As this is again 
a linear stochastic differential equation, we can follow 
the same strategy as in the present work and investigate 
the response to a temperature quench by analyzing exact 
expressions. On the other hand, we can also extend our 
study to systems that are of direct relevance for thin 
film growth: the Kardar-Parisi-Zhan g (K PZ) [§] and the 
conserved KPZ universality classes (lOL |Tl|. As exact 
expressions for these non-linear systems are not available, 
we plan to integrate these equations numerically and to 
simulate quenches in microscopic models belonging to the 
same universality classes. 
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Appendix 

Here, we show how to compute the first crossover time 
ti in the limit of infinite L (see also ^3^). Since this point 
is defined as the intersection of the RD regime ( W"^ = Dt) 
and the EW regime {W"^ = At^^^), we have 

h = (A/D)^ (16) 

so that the problem reduces to finding the amplitude A 
associated with EW growth. For L ^ oo , the sum in 



expression ([3]) can be replaced by an integral, which can 
be computed to extract A. Further simplification occurs 
if we consider dtW^ = D /^^ e^^"*^ d9/2TT instead. Im- 
posing the ansatz = At^^^, we arrive at 

r e-2-*e^d0. (17) 

7^ J-7T 

Transforming to ^ = \/2iA9, we have 

A^-=_/ e-^'dt (18) 

So, as t — > oo (or, to be precise, L ^ t ^ we arrive 

at = ^ ^gij 

\D J 2'Kv V 

with T = l/27r. 

Of course, this approach can also be used to extract 
t-i. Equating this W"^ (i.e., At^l"^) to the saturation W"^ 
(i.e., {DLj^-K^v) J2n'^'^ = DL/2Av), we arrive at 



This provides 




- = log ^7^1/12. 

We should emphasize again that these results are exact 
in the L ^ oo limit. 
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